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ABSTRACT 

A 227 ksec Chandra Observatory X-ray image of the hot plasma in the Hydra A cluster has revealed an 
extensive cavity system. The system was created by a continuous outflow or a series of bursts from the nucleus 
of the central galaxy over the past 200-500 Myr. The cavities have displaced 10% of the plasma within a 300 
kpc radius of the central galaxy, creating a swiss-cheese-like topology in the hot gas. The surface brightness 
decrements are consistent with empty cavities oriented within 40 degrees of the plane of the sky. The outflow 
has deposited upward of 10 61 erg into the cluster gas, most of which was propelled beyond the inner ~ 100 kpc 
cooling region. The supermassive black hole has accreted at a rate of approximately 0.1-0.25 M Q yr" 1 over 
this time frame, which is a small fraction of the Eddington rate of a ~ 10 9 M black hole, but is dramatically 
larger than the Bondi rate. Given the previous evidence for a circumnuclear disk of cold gas in Hydra A, 
these results are consistent with the AGN being powered primarily by infalling cold gas. The cavity system is 
shadowed perfectly by 330 MHz radio emission. Such low frequency synchrotron emission may be an excellent 
proxy for X-ray cavities and thus the total energy liberated by the supermassive black hole. 
Subject headings: cooling flows — galaxies: clusters: general — galaxies: elliptical and lenticular, cD — 
intergalactic medium — X-rays: galaxies 



1. INTRODUCTION 

The realization that powerful outbursts from central active 
galactic nuclei (AGN) can release upward of 10 61 erg into the 
intracluster medium (ICM) points to a common solution to 
several heati ng problems associated with clusters and group s 
of galaxies dMcNamara et all 120051: iNulsen et all l2005albh . 
AGN feedback has become the most promising mechanism 
for regulating the cooling of the hot gas in galaxies and clus- 
ters, and thus limiting the growth of ellipti cal galaxies and 
their attendant supermas sive black holes dVoit & Donahue! 
120051: iBenson et alj|2003l) . Furthermore, energy released by 
AGN could be a major source of excess en tropy ("preheat- 
ing") in th e hot gas in groups and clusters dWu et all 120001: 
IVoitl 120041) . and AGN outbursts may thread clusters with 
large-scale magnetic fields. 

Cavities and shock fronts provide a reliable lower limit to 
the energy output of AGN ind ependent of the magnetic field 
and particle content of the jet (Bi rzan et al.ll2004l iDunn et all 
l2005t iDeY oung 2006). The existence of well-defined cavi- 
ties in the hot ICM around radio lobes implies that the mag- 
netic field and relativistic particles are largely confined to the 
volume of the cavity such that most of the energy output 
of the AGN is included (radiation losses are usually negli- 
gible in these systems). This energy provides a strict lower 
limit to the gravitational binding energy released by accre- 
tion onto the embedded supermassive black hole, and hence 
the mean growth rate of the hole over the life of the outburst 
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dMcNamara etal.1 12005: Nuls en et al.1 l2005at iRaffertv et all 
2006). 

This geometric approach is limited at present by the Chan- 
dra X-ray Observatory's ability to resolve cavity boundaries 
required to measure their sizes (volumes) and the surround- 
ing pressures. This provides a largely parameter free esti- 
mate of the work done by radio jets as they inflate the cav- 
ities against the surrounding gas pressure (pV). The method 
assumes the cavities are close to pressure balance with the 
surrounding gas, which is justified by the absence of bound- 
ing strong shocks, their high detection rate and advanced ages 
(Birzan et al. 2004; iDunn et alj[2005b . The total free energy 
(enthalpy) per cavity is expected to lie in the range 2.5 pV - 
4pV if the pressure is dominated by gas, or 2pV if it is domi- 
nated by magnetic field. Thus, the free energy depends on the 
equation of state of the (radio) plasma filling the cavities. The 
principal sources of uncertainty in the enthalpy estimates are 
the volume measurements (geometry) and the as yet unknown 
makeup (equation of state) of the plasma inside the cavities, 
which combined amount to uncertainties of factors of a few. 

The total energy pumped into the ICM by AGN is par- 
titioned into radio emission (which is generally negligible), 
cavity energy (3pV for relativistic gas), and "shock energy." 
Here "shock energy" refers to the total work done by the 
expanding cavities on the surrounding gas, which appears 
as excess thermal, kinetic, and potential energy in the ICM. 
Its effect is most evident where the expanding cavities drive 
a detectable shock. Measurements of these quantities have 
shown that the internal cavity pressures generally exceed the 
radio equipartition pressures by more than an order of mag- 
nitude and that jet (mechanical) power often exceeds syn- 
chrotron power by several orders of magnitude. The unavoid- 
able conclusion then is that radio jets possess much larger 
mechanical powers than expected, with potentially dramatic 
con sequences for qu enching cooling flows and cluster heat- 
ing dPedlaret alj|1990h . 

Whether AGN heating of the ICM occurs violently through 
shocks or outflows, or gently through bubbling, vorticity, etc., 
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FIG . 1 . — Chandra 0.5-7.0 keV mosaic of the Hydra A field. The image has been background-subtracted and exposure-corrected before smoothing with a 1 .5" 
Gaussian. The horizontal bar shows the spatial scale in the image. A logarithmic color table has been used to highlight low surface brightness features. 



or both is poorly understood. Most known cavity systems 
are associated with Fanaroff-Riley class I radio sources, for 
which jet ram pressure is probably negligible. In this instance 
the pr essure within the X-ray cavities is expected to be uni- 
form dHeinz et al.ll 19981) . and thus the cavity pressure is well- 
defined by the external pressure. The ratio of cavity enthalpy 
to shock energy is governed chiefly by the history of power in- 
put to the cavity. Therefore, the nature and history of energy 
injection can be evaluated using the cavity properties and the 
surrounding shock front. Here we examine this and other is- 
sues, using an analysis of the enormous system of X-ray cav- 
ities embedded in a "cocoon" shock in the Hydra A cluster 
(McN amara et al.ll2000tlDavid et al.l20"0ltlNulsen et al.l2002l 
2005b), and we discuss its consequences for cluster heating 
by the central radio source, and for accretion onto the central 
supermassive black hole. 

In the succeeding analysis, we have adopted a redshift for 
Hydra A of z = 0.0538 and a flat ACDM cosmology with 
H()=70 km s" 1 Mpc -1 , $7^=0. 3, and fl\=0.7. These assump- 



tions yield a luminosity distance of D = 240 Mpc, an angular 
diameter distance of 216 Mpc, and a linear scale of 1.05 kpc 
per arcsec. 

2. OBSERVATIONS AND DATA REDUCTION 

Hydra A has been observed four times by Chandra with 
the ACIS detector in imaging mode for a total exposure of 
240 ksec. Three of the observations were made with ACIS-S 
detector at the aimpoint (ObsIDs 576, 4969, 4970) and one 
(ObsID 575) using the ACIS-I detector. The two shorter ex- 
posures were taken early in the Chandra mission when the 
ACIS detector was operated at the higher focal plane temper- 
ature of -110° C. Each dataset was individually reprocessed 
using CIAO 3.2 and calibration files available in CALDB 3.0. 
Standard screening was applied to the event files to remove 
bad grades and pixels. Two of the datasets (ObsIDs 575 and 
4969) exhibited strong background flares and required addi- 
tional filtering. The final, combined exposure time for all four 
datasets was 227 ksec. 



Supercavities in Hydra A 




FIG. 2. — High-contrast residual map after subtraction of a smooth, elliptical beta model fit to the cluster surface brightness profile. The positions of the three 
pairs of cavities discussed in the text are indicated. The position of the central radio source is indicated by the light blue dot. 



After cleaning, the event files for ObsID's 575, 576, and 
4969 were reprojected to a common tangent point on the sky 
with ObsID 4970. A number of bright point sources in Hydra 
A were used to test for potential astrometry errors in the repro- 
jected sky coordinates between fields. Individual background 
event files were created for each ObsID from the standard 
A CIS blank-sky event fil es following the procedure described 
in lVikhlinin et al.1 (120051) and reprojected to the common tan- 
gent for ObsID 4970. Using the reprojected event files, source 
and background counts images were created for each CCD in 
each observation using standard CIAO tools. Matching expo- 
sure maps were created for each chip using the CIAO tool 
expmap. The resulting images for individual CCDs were 
then added to form mosaics combining all four datasets. The 
resulting background-subtracted, exposure-corrected mosaic 
is shown in Figure Q]centered on the radio source 3C 218. All 
detected point sources were excised prior to extracting spectra 
and surface brightness profiles. 

Due to the higher quality calibration at -120° C, spectral 
analysis was restricted to the two longer, more recent expo- 
sures (4969 and 4979). Collectively, ObsIDs 575 and 576 
comprise 17% of the total exposure and neglecting them has 
little or no effect on the final spectral fits. Count weighted 
detector response (RMFs) and effective area (ARFs) files 
were created for extraction regions using the CIAO tools 
mkacisrmf and mkwarf, including the temporal, spectral, 
and spatial dependences of the ACIS filter contaminant. 



3. CLUSTER-SCALE RADIO CAVITIES & FILLING FACTOR 

In order to highlight faint cavity structures, an elliptical beta 
model was fit to the azimuthally averaged surface brightness 
and subtracted. The residual image shown in Figure [2] re- 
veals at least three pairs of X-ray cavities emanating from the 
central cD galaxy distributed with a "swiss cheese" topology 
throughout the inner region of the cluster. 

The inner pair of cavities (label ed A and B in F i gure [3 
are well-known se t described b y McNamara et al. (2000); 
iDavid et al.1 (1200 ll) : iNulsen et all (120021) . These cavities are 
each 40-50 kpc in diameter and centered roughly 25 kpc 
north-south of the central AGN. We also note that the en- 
hanced X-ray emission, relative to the underlying smooth 
model, immediately to the west of cavit y A corresponds to t he 
region of cooler gas seen previously in Nul sen et al.l (120021) . 

A second set of "ghost" cavities (labeled C and D) can be 
seen beyond the inner pair. The north ern-most correspond s 
to the cavity "C" feature discussed in Nuls en et al.l (|2005b). 
This asymmetrical pair of cavities has diameters of 50 kpc 
and 80 kpc, projected 100 kpc and 60 kpc distant from the 
north and south cavities, respectively. Finally, two very large 
cavities (E and F) are projected at radii of 225 kpc and 100 
kpc to the north and south, respectively. These cavities have 
diameters of 120-200 kpc and fall just inside a sharp edge in 
the surface brightness w hich can also be seen in Figure|2] This 
edge was interpreted bv lNulsen et al . (2005b) as a shock front 
associated with a 140 Myr old outburst. The estimated Mach 
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TABLE 1 
Cavity Properties in Hydra A 





a a 


b b 


R' 


P V 
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P h 


Cavity 


(kpc) 


(kpc) 


(kpc) 


(10 58 ergs) 


(10 7 yr) 


(10 7 yr) 


(10 7 yr) 


(10 7 yr) 


(10 44 ergs s-') 



A 20.5 12.4 24.9 8.4 

B 21.0 12.3 25.6 8.2 

C 47.2 31.5 100.8 29.3 

D 29.0 20.9 59.3 14.8 

E 105.0 99.7 225.6 247.9 

F 67.7 50.1 104.3 102.3 



2.7 4.3 8.6 5.2 2.1 

2.7 4.4 8.7 5.3 2.0 

10.2 23.0 26.8 20.0 1.9 
6.1 13.3 16.4 11.9 1.6 

22.3 51.7 65.3 46.4 6.8 
10.5 20.2 33.5 21.4 6.1 



a Projected semimajor axis of the cavity. b Projected semiminor axis of the cavity. c Projected distance from the cavity center to the radio core. d Estimated cavity 
age based on the sound speed, c s , in the gas. e Estimated cavity age assuming the cavity rises buoyantly.' Estimated cavity age based on the refill timescale. g Mean 
cavity age for all three estimates. 11 Cavity power is calculated assuming 4pV of energy per cavity and the mean timescale for the age of the cavity (f). 



number for this shock is roughly 1.2, although the data are 
not of sufficient quality to detect the 20% jump in temperature 
expected for such a shock. 

This enormous cavity system fills at least ~10% of the clus- 
ter's volume within a 300 kpc radius (see below). Table[T]de- 
scribes the cavities in the system. All appear to be connected, 
raising the possibilities that they are merging bubbles gener- 
ated by a series of outbursts, or several attached cavities that 
have been continuously powered by the AGN over the past 
several hundred million years. 

4. CAVITY GEOMETRY & ENERGETICS 

The deprojected gas densities and temperatures surround- 
ing the c avities were determined using the t echniques de- 
scribed in IWise. Mc Namara. & Mu rray! (120041) and are con- 
sistent with those of iDavid et al.l d2001l) . Values for the sur- 
rounding pressures were combined with the volume estimates 
to calculate the IpV energies of each cavity listed in TableQ] 
Summing the pV values for all six cavities from Table Q] the 
total work done by the jet is 4.1 x 10 60 ergs, correspond- 
ing to a total free energy or enthalpy {2pV - 4-pV) of 0.8 - 
1 .6 x 10 61 ergs. This range is comparable to the energetic out- 
put of other supercavity sys tems MS0735. 6+7421: 8 x 10 61 
erg (McNamara et al.l 120051) . and Hercules A: 3 x 10 61 erg 
(Nulse net alJl2005alK It is also comparable to the ~ 1 x 10 61 
erg lNulsen et al (2005b) found for the energy associated with 
the large-scale shock in Hydra A. 

Inspection of the cavity energies in Table Q] shows that the 
outer, larger cavities contain more than an order of magnitude 
more energy than the original A and B cavity systems. Most 
of this energy is being deposit ed into the ICM outside the 100 
kpc c ooling radius in Hydra A (IDavid et al.l2001l:lBrrzan et al.1 
2004). Less than 10% of the total free energy associated with 
the cavity system resides inside the cluster's cooling radius. 

A number of factors contribute to the uncertainty in these 
energy estimates including geometric uncertainties in the size 
and shape of the cavities, the "filling factor" or fraction of 
the volume from which X-ray emitting gas is excluded, and 
the locations of the cavities along the line-of-sight through 
the cluster. Since these factors also affect the observed X- 
ray surface brightness profile over the cavities, these profiles 
can be used to constrain these uncertainties and to test our as- 
sumptions about the cavities. For example, the presence of an 
empty cavity within the cluster emitting volume will produce 
a decrement in the X-ray surface brightness that depends on 
the size of the cavity and its depth in the cluster. By mea- 
suring these cavity decrements, we can therefore set limits on 
the placement of the cavity along our line-of-sight through the 



cluster. 

Table|2]gives the decrement for each cavity, y, defined to be 
the ratio of the surface brightness measured at the center of 
each cavity to the surface brightness at the same radius in the 
"undisturbed" ICM. The undisturbed ISM is modeled using 
the the double beta model fit to the azimuthally av eraged sur- 
face br ightness profile (excluding the cavities) of David et al. 
d2001l) . This model is described by values of /3i = 0.686, 
02 = 0.907, core radii of 27.7 kpc and 235.6 kpc, respectively, 
and a normalization for the more extended, second beta model 
component relative to the first of 0.0418. To estimate the er- 
ror in the decrement, <5y, we determined the distribution of y 
values in circular apertures 5 arcsec in radius, centered on the 
cavity centers. The decrement errors, 8y, quoted in Table [2] 
represent the la widths of these distributions. 

If we assume that the width of the cavity along the line of 
sight is 2r e ff = 2(ab) l l 2 and the underlying surface brightness 
profile is given by the double beta model fit described above, 
we can determine the distance, z, of the cavity along the line- 
of-sight. This distance is measured from a plane perpendicu- 
lar to the line-of-sight which passes through the central AGN. 
Table |2]lists the values of z for each cavity required to repro- 
duce the observed surface brightness decrement y. One sigma 
ranges for z, z m in - Zmax are also given. Missing entries fall 
outside the range of surface brightness allowed by this model. 
These values are insensitive to details o f the surface brightn ess 
model. Using the single beta model of Davi d et all (1 19901) for 
the outer cavities, or a single beta model for the Chandra data 
at small radii for the inner cavities gives very similar results. 

With the possible exception of cavity A, the results in Ta- 
ble |2] show that the the cavities must be almost devoid of X- 
ray emitting gas at the ambient temperature and density and 
lie close to the same plane as the AGN. If their filling factors 
are less than unity, then they would need to lie even closer to 
the plane of the AGN. Thus, the surface brightness deficit for 
most of the cavities is consistent with our key assumptions. 
Only the deficit of 0. 1 for cavity A is marginally smaller than 
expected. Its elongation along the northern radio "jet" sug- 
gests that cavity A may be prolate; however, replacing the ef- 
fective radius with the semiminor axis in the model only shifts 
it a little closer to the plane of the AGN, from 45 kpc to 39 
kpc. At a distance of 45 kpc, it makes an angle of ~ 60° to 
the plane of the sky, well within reasonable bounds. However, 
the alignment of the inner radio lobes (cavities A and B) in the 
plane of the sky suggests that they should make the same an- 
gle to the plane of the sky. Since the other cavities all make 
smaller angles to the plane of the sky, cavity A appears to be 
the anomaly. Its low decrement may be due to dense, cool gas 
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TABLE 2 
Cavity Decrements in Hydra A 







R 


y" 


8f 


z d 


£min 


-'■max 


atan(z/i?) 


Cavity 


(kpc) 


(kpc) 






(kpc) 


(kpc) 


(kpc) 


(deg) 



A 15.97 24.9 0.90 0.06 45 35 65 61 

B 16.06 25.5 0.67 0.04 19 15 22 37 

C 38.58 100.8 0.71 0.08 26 - 70 14 

D 24.62 59.3 0.59 0.05 - 

E 102.34 225.6 0.62 0.15 101 - 193 24 

F 58.26 104.3 0.71 0.09 80 43 118 37 



a For ellipsoidal cavities, the effective radius is taken to be r c ff = (afo) 1 / 2 , with values taken from TablefT] b The cavity decrement as discussed in the text, expressed 
as the ratio of the surface brightness at the center of the cavity relative to the underlying smooth surface brightness model, y = Scav/Ssmo- Values less than one 
show decrements. Estimate of the lcr error in the cavity decrement based on the width of the distribution of y values as measured in a 5 arcsec radius, circular 
aperture centered on the cavity center. d Distance from the plane of the sky containing the AGN along the line-of-sight to the cluster. 



lying over the cavity. In any case, cavity A contributes little 
to the total energy. 

By contrast, the deficit for cavity D is too large to be ex- 
plained by our model. Cavity D lie s where the south ern radio 
lobe appears to fold back on itself (lLane et al.ll2004l also Fig- 
ure |3), so that our line-of-sight may be nearly tangent to the 
lobe axis there. In that case, it would be reasonable to expect 
its central depth to exceed 2r e ff. If true, the pV value quoted 
for cavity D in Table Q~]may be an underestimate. 

5. CAVITY AGES AND POWERS 

Analogous to the well-known 1 .4 Ghz radio emission filling 
the innermost cavities, low frequency radio observations of 
Hydra A show that these large-scale cavities also contain rel- 
ativistic plasma strongly suggesting a connection to the AGN 
located in the nucleus of the central cD galaxy. This point is 
illustr ated in Figu re [3] which overlays the 330 Mhz VLA im- 
age of Lan e et all (120041) with the residual X-ray image shown 
in Figure[2] The existence of multiple cavities implies that the 
central AGN is periodically rejuvenated producing a series of 
outbursts. 

We can gauge the timescale over which this activity has oc- 
curred by deter mining the ages of th ese cavities. Following 
the discussion in Bir zan et al.l (120041) . we have estimated the 
age of each cavity in three ways. First, the age of the cavity 
was estimated as the time required for it to reach its projected 
location assuming it traveled at the sound speed. Alterna- 
tively, the age was estimated as the time for the cavity to rise 
buoyantly to its present l ocation. Finally, the age was taken 
to be the refill timescale (McNamara et al. 2000; Nulsen et al. 
2002). Each of these timescales is tabulated in Table Q] along 
with the mean cavity age, (f), for the three methods. Cav- 
ity ages based on the sound speed typically yield the small- 
est timescales whi le refill timescales tend to be the largest (cf 
iBirzan et al.ll2004|) . We find m ean ages of ~ 50 M yr for cav- 
ities A and B, consistent with (Birzan et al. 2004) and values 
ranging between ~ 100 - 200 Myr for cavities C and D. Cav- 
ity E could be up to almost 500 Myr old. 

Combining the values of pV given in Table Q] with these 
age estimates, we can calculate the mean instantaneous power 
of each cavity as P cav = 4pV / (t). Here, we have explicitly 
assumed the cavity is filled predominantly with relativistic 
plasma, making its enthalpy 4pV. The resulting values for 
f C av are presented in TableQ] Equating P cav with the output of 
the AGN, the data for the four inner cavities imply that a pair 
of cavities is created in Hydra A every ~ 50 - 100 Myr for a 
mean AGN power of 2 x 10 44 ergs s -1 . On the face of it, the 
larger, outermost cavities would seem to be associated with 



an earlier, more energetic state of activity which occurred 200 

- 500 Myr ago with an associated AGN power output of ~ 6 

- 7 x 10 44 ergs s" 1 . 

6. SHOCK FRONT AGE AND POWER 

By fitting simple shock models to the surfa ce brightness 
disconti nuity in Hydra A observed at ~ 210 kpc. [Nulsen et al.1 
(2005b) estimated an age for the shock of 140 Myr, smaller 
than the shortest age of 220 Myr for cavity E found here. At 
its closest point to the AGN, the projected distance from the 
shock front to the AGN should be very nearly equal to the ac- 
tual distance, so that the age estimate based on the shock is 
likely to be the most accurate. Given their relative locations, 
it is not surprising that an age based on the closest points of 
the (supersonic) shock front is smaller than any estimate in 
Table Q] for the age of cavity E. This result is a reminder that 
a momentum dominated jet can excavate a remote cavity in 
considerably less time than it would take for the cavity to rise 
buoyantly from the AGN. It also illustrates that cavity based 
ages may be too high, and the total power may be underesti- 
mated. 

The preceding interpretation presupposes that the cavities 
do not share energy, whereas the radio and X-ray maps show 
that they might be interlinked (Figure |3). Furthermore, if the 
cavities are disconnected, then the jet cannot supply power 
to the outermost cavities that are chiefly responsible for driv- 
ing the shock. We should then expect the shock to separate 
rapidly from those cavities. For example, cavity E needs to 
be overpressured by a factor cl ose to 2 in order to drive the 
Mach 1.34 shock to the north dNulsen et al.ll2005bl) . If it is 
dominated by relativistic gas, then it can only expand adia- 
batically by about 19 kpc from its current radius, ~ 100 kpc, 
before coming to local pressure equilibrium. At a speed of 
356 km s" 1 (see below), that would take ~ 50 Myr, leaving a 
narrow window after it is disconnected from the jet in which 
cavity E can continue to drive the shock. A more plausible 
alternative is that the cavities remain interconnected. In that 
case, energy in the inner cavities is only in transit and cannot 
be used to estimate jet power (in the absence of a much more 
sophisticated model). Only the outermost cavities, which are 
the destination of the energy carried by the jet, can be used 
to estimate its power. If the northern jet remains connected 
to cavity E, it is also easier to explain how the cavity remains 
close to the shock front. 

In these circumstances, the proximity of the shock to cav- 
ity E suggests another estimate for the power of the north- 
ern jet. Radio losses are negligible, so that jet power is di- 
vided between cavity thermal energy and p dV work done by 
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FIG. 3. — Composite color image of Hydra A which illustrates the close connection between the observed, large X-ray cavity system (shown in blue) and the 
low frequency, 330 Mhz radio emission (shown in green). The X-ray emission corresponds to the residual image shown in Figureff] The 330 Mhz radio data is 
from lLane etaT] flOOl) . The familiar 1 .4 GHz VLA image of Hydra A is also shown in the core in yellow. 



the expanding cavity, i.e., Pj et = dE cav /dt + pdV cav /dt. As- 
suming that evolution of the cavity and shock is approxi- 
mately self-similar, the radius of the cavity scales with the 
radius of the shock, giving dR CRV /dt ~ v S hockCRcav/^shock) — 
356 km s" 1 , forfl cav = 1 00 kpc, fl sho c k = 35 kpc, kT = 3.2 keV 
for the unshocke d gas dDavid et al.ll2001l) and a shock Mach 
number of 1.34 dNulsen et all l2005bl) . Thus pdV^/dt ~ 
47r/?/?;? av v s hock/fishock — 6.4 x 10 44 erg s~\ for the parameters 
above and a cavity pressure of 1.5 x 10 -11 erg cm -3 . If the 
pressure of the cavity remains constant, including the increase 
in thermal energy boosts this by 7/(7- 1), where 7 is the 
ratio of specific heats for the lobe plasma. Allowing for a 
decrease in cavity pressure, p oc K^L, modifies the boost to 
7/(7- l)-77/3. For self-similar growth, the cavity pressure 
follows the external pressure and 77 ~ 2 from above, so that, 
for 7 = 4/3 in the cavity and we get fj et ~ 2 x 10 45 erg s -1 . Al- 



though evolution of the cavity and shock is unlikely to be ex- 
actly self-similar, this estimate is better determined than those 
relying on cavity rise times. 

7. HISTORY OF THE OUTBURST 

Jet power is divided between thermal energy in the cavities, 
~ pV cav , and the work done by the cavities on the surrounding 
gas as they expand, i.e., J pdV cdv over the history of the cavi- 
ties (radiation losses are negligible). The ratio pV cav / j pdV cdv 
is then a diagnostic for the history of an AGN outburst. In- 
creasing cavity volume and decreasing external pressure both 
tend to make the cavity pressure decrease with time, so that 
jet power would need to increase rapidly to maintain constant 
pressure. Thus, unless the jet power is increasing rapidly, 
pKav/ / pdV c . dv < 1. This ratio is smallest when the jet power 
is declining. 
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The observations do not give a clear indication of whether 
the jets are steady, ramping up, or declining. From § @] the 
total of pVcav for the cavities is ~ 4 x 10 60 erg, which can be 
compared to the t otal energy of ~ 9 x 10 60 erg required for 
the shock model dNulsen et al.ll2005bh . In principle, the lat- 
ter includes both the thermal energy of the cavities and the 
work done by them. This shock model is spherical, is fitted 
to the part of the shock front closest to the AGN, and assumes 
a single, explosive energy injection event at the cluster cen- 
ter, rather than continuous energy injection via jets. These 
approximations tend to underestimate the thermal energy of 
the cavities. Nevertheless, the numbers suggest that the ratio 
pVcav/ J pdV cav is not much smaller than unity. On this ba- 
sis, the jet power does not appear to be declining rapidly at 
present. 

On the other hand, the sharp bend in the northern jet at cav- 
ity C indicate s that the jet may no t be momentum dominated 
at larger radii dNulsen et al.ll2005bl) . If true, the long rise time 
of cavity E compared to the age of the shock would require 
that it was partially inflated in an earlier outburst. The ear- 
lier outburst contributes to pV caw , but not the current shock, so 
that pVcav/ / pdVcw would be larger than expected for a sin- 
gle outburst. Jet power may then be declining more rapidly 
than suggested in the preceding paragraph. We note that the 
shape of the radio source and the properties of cavity D are 
also consistent with a sharp bend in the southern jet there (see 

There are other indications that the situation is more com- 
plex than we might have imagined. For example, cavities C 
and D may be kept open by the forces required to deflect the 
jet. If the total jet power is ~ 2 x 10 45 erg s" 1 (§ [5]l and half 
of this passes through the northern jet, then its minimum mo- 
mentum flux is F = Pj/c ~ 3 x 10 34 dyne. The force required 
to deflect the jet through an angle 6 is then 2F sin 6/2 ~ F, for 
a deflection of 60°. If the outer wall of cavity C is cylindrical, 
with radius of curvature r ~ 50 kpc and depth 2r (cf. Table[T), 
the minimum force per unit area on the cavity wall needed to 
deflect the jet is P def ~ F/(nr 2 /3) ~ 1.3 x 10" 12 erg cm -3 . 
This value is increased by a factor of 2c/ v, if the power of the 
jet is primarily kinetic energy and its bulk speed is v <C c. The 
gas pressure around cavity C is ~ 5 x 10" 11 erg cm -3 . Thus, 
the force required to deflect the jet could keep cavity C open 
if v~ 0.1c. 

Another issue for the history of the outburst is formation 
of the innermost cavities. In order to open them and/or keep 
them open, the pressure in cavities A and B must equal or ex- 
ceed the pressure of the local ICM. Furthermore, any recent 
pressure increase (change) must be slow in order to avoid cre- 
ating new shocks around these cavities. Thus, if the inner- 
most cavities were created by a recent boost in AGN power, 
then it needs to have been gradual, over a timescale signif- 
icantly longer than their sound crossing time of ~ 30 Myr 
(Table Q]). The pressure in the ICM around cavities A and B 
is also an order of magnitude greater than that around the out- 
ermost cavities. Since the pressure changes are gradual, we 
can apply Bernoulli's theorem to the jet (whether relativistic 



or not) and the large pressure change requires the bulk speed 
of the jet flow to become highly supersonic (until a shock is 
encountered). 

In summary, it is unclear whether the existence of multiple 
cavities requires multiple outbursts. The issues raised here 
highlight the need for more detailed models in order to inves- 
tigate the history of this system. 

8. DISCUSSION 

A deep Chandra image of Hydra A has revealed a complex 
cavity system in the central 300 kpc of the Hydra A cluster in- 
dicating a long and equally complex AGN outburst history en- 
ergetically equivalent to that of a powerful quasar. The AGN 
has been active either continuously or intermittently for the 
past several hundred million years, as it has deposited several 
10 61 erg into the surrounding cluster gas. The corresponding 
jet power 2 x 10 45 erg s" 1 is enough energy to offset radia- 
tion losses a ssociated with a se veral hundred M yr" 1 cool- 
ing flow (eg. iDavid et alj.l2001h . and to heat the gas well be- 
yond the cooling region. The cavities occupy ~ 10 percent of 
the volume within 300 kpc of the center, and they lie within 
about 40 degrees of the plane of the sky. The cavities are filled 
entirely by 330 MHz radio emission, demonstrating that low 
frequency radio emission faithfully traces the energy of the 
outburst. 

The association between the observed radio emission and 
the X-ray cavities implies that these cavities were created 
and powered by accretion onto a supermassive black hole 
embedded in the cD. The inferred mass of the black hole 
based on the relationships between black hole mass, stellar 
velocity dispersi on, and cD galaxy l uminosity corresponds to 
- 9 x 10 s M Q jRaffertvet alJl2006h . The power, P, liber- 
ated over the lifetime of the event implies an accretion rate of 
P/ec 2 ~ 0.05-0.25 M© yr" 1 , where we have adopted e = 0.1 
as the mass to energy conversion efficiency. The range of 
values depends on the adopted age and total energy of the 
outburst. By comparison, the Eddington rate for a black 
hole of this mass is ~ 20 M yr" 1 , implying the accretion 
driving the outburst is in the strongly sub-Eddington regime. 
The Bondi accretion rate implie d by the black hole m ass and 
central gas density taken from ( Raffer tv et al J 120061) is only 
~4x 10" 4 M Q yr" 1 , which is too small to power an outburst 
of this magnitude. On the other hand, there is no shortage of 
cold fuel near the nucleus. Approximately 10 4 - 10 7 M of 
neutral hydrogen is present in a circ umnuclear di sk or torus 
located within 5-30 pc of the nucleus dTaylorll 19961) . This ma- 
terial is enough to fuel the outburst for several tens of thou- 
sands to several hundred million years. Additional hydrogen 
and CO is likely pre sent in a larger cir cumnuclear disk of gas 
and star formation ( McNamaral l 1 9951) in amounts that do not 
exceed - 10 9 M (Edge 2001). 
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